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Many outbreaks of A(HlNl)pdm09 influenza have occurred in schools with a high population density. 
Containment of school outbreaks is predicted to help mitigate pandemic influenza. Understanding disease 
transmission characteristics within the school setting is critical to implementing effective control measures. 
Based on a school outbreak survey, we found almost all (93.7%) disease transmission occurred within a 
single grade, only 6.3% crossed grades. Transmissions originating from freshmen exhibited a star-shaped 
network; other grades exhibited branch- or line-shaped networks, indicating freshmen have higher activity 
and are more likely to cause infection. Rq for freshmen, calculated as 2.04, estimated as 2.76, was greater than 
for other grades (P < 0.01). Without intervention, the estimated number of cases was much greater when the 
outbreak was initiated by freshmen than by other grades. Furthermore, the estimated number of cases 
required to be under quarantine and isolation for freshmen was less than that of equivalent other grades. So 
we concluded that different grades have different transmission mode. Freshmen were the main facilitators of 
the spread of A(HlNl)pdm09 influenza during this school outbreak, so control measures (e.g. close contact 
isolation) priority used for freshmen would likely have effectively reduced spread of influenza in school 
settings. 



The 2009 HlNl influenza virus pandemic (A(HlNl)pdm09) emerged in Mexico and rapidly spread globally, 
resulting in millions of cases and approximately 18,000 deaths in nearly 200 countries^^'^^. Surveillance data 
from many countries have shown that school- aged children and young adults most likely to be infected by 
A(HlNl)pdm09 influenza^^'^^'^^'^^. Many outbreaks occurred in schools, presumably because of a higher density 
of susceptible persons^°'^^'^^'^°'^^"^\ Surveillance data on A(HlNl)pdm09 influenza in mainland China from May 7 
to November 30, 2009 indicates that 39% (27,806/71,665) of reported A(HlNl)pdm09 cases were reported to be 
part of a recognized outbreak and 94% of those cases associated with an outbreak were 5-24 years of age^^. Based 
on surveillance data from Beijing city, 81% (47/58) of outbreaks with 10 or more laboratory- confirmed cases 
occurred in schools, and most of school outbreak was initiated by freshmen. Studies have also suggested that 
school-aged individuals facilitated the introduction and spread of A(HlNl)pdm09 influenza in households''^. 
Curbing school outbreaks would therefore help to mitigate pandemic influenza. Understanding disease trans- 
mission characteristics within school settings would therefore be helpful in defining effective control measures for 
A(HlNl)pdm09 influenza. 

Compared to Western countries, there have been few descriptions of the transmission dynamics of 
A(HlNl)pdm09 influenza in school outbreaks in mainland China, where campus quarantines (fengxiao) rather 
than school closures represent the primary public health intervention^'^'^°'^^'^^'^^''°'^\ Studies have suggested that 
social networks and population structures affect the spread of communicable diseases and the selection of optimal 
control strategies^'^'^^''\ To identify optimal control interventions, we used an epidemiological survey and disease 
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transmission dynamic model to explore the disease transmission 
characteristics among students in a school involved in a outbreak 
of A(HlNl)pdm09 influenza. 

Results 

Epidemiological Features. The epidemiological investigation 
showed that this outbreak was initiated by one infected freshman. 
A total of 226 confirmed cases were found, comprising 212 students 
and 14 teachers. Among confirmed cases, 198 were male (93.3%) and 
14 were female (6.7%). This study focused on disease transmission 
among the student population, so 212 students data were used to 
analysis below. Of the 212 infected students, 24 (11%) were also 
excluded from the transmission chain analysis because they could 
not recall the source of their infection or contacts. Among 212 
infected students, 118 (63%) were freshmen, 65 (35%) sopho- 
mores, 16 (9%) juniors, and 13 (7%) seniors (Figure 1). The infec- 
tion rate was 6.8% among freshmen, 3.0% among sophomores, 0.9% 
among juniors and 0.95% seniors. Freshmen were the most affected 
group, with the highest infection rate among the four grades (Chi 
square test, p < 0.0001). 

Transmission Characteristics. The transmission tree of A(HINI) 
pdm09 influenza among the confirmed cases is shown in Figure 2a. A 
total of 142 transmission chains were identified, comprising 133 
(93.7%) within-grade chains and 9 (6.3%) chains across grades 
(Figure 2b). Moreover, disease transmission with a freshman as the 
infection source exhibited star- shaped networks; and transmission 
with students from the other three grades as the infection source 
exhibited branch- or line-shaped networks (Figure 2a). This 
indicates that freshmen were the important determinant of social 
mixing between the students. An investigation of whether fresh- 
men affected transmission dynamics of A(HlNl)pdm09 influenza 
in the school was performed by calculating basic reproductive 
numbers (Rq) from the different grades. The basic reproductive 
number of an infection, Rq, is defined as the expected number of 
secondary infectious cases generated by an average infectious case in 
an entirely susceptible population^^, and it was calculated from 
epidemiological data using the average number of individuals 
directly infected by a confirmed case during her or his entire 
infectious period'*^. The Rq values were 2.04, 1.18, 1.40, and 1.00 
for freshmen, sophomores, juniors and seniors, respectively. The 
95% confidence intervals of RO are [1.61, 2.45], [0.99, 1.38], 
[0.7199, 2.08] and [1, 1] for freshmen, sophomores, juniors and 
seniors, respectively. Because the numbers of secondary infectious 
cases resulting from seniors were all 1 based on our epidemiological 
investigation, so the confidence interval of the RO for seniors is [1, 
l].The RO value of the seniors was calculated based on data from the 
epidemiological investigation, the standard estimate of Rq is always 
an underestimate of the true value^°. So the actual Rq of seniors 



should be higher than 1.0. The Rq value for freshmen was signifi- 
cantly higher than those for the other grades (single factor analysis of 
variance, P < 0.01; Figure 3). 

Estimated Transmission Dynamics. Initial dynamics of disease 
transmission among the different grades were reproduced using 
the constructed model. The parameters were estimated using our 
model: oli = 1.1, as = 0.8, = 0.70, 0^4 = 0.69, as = 0.30 and C 
= 0.05. The estimated Rq values determined by the model were 2.76, 
2.01, 1.76, and 1.73 for freshmen, sophomores, juniors and seniors, 
respectively. The 95% confidence intervals of RO are [2.73, 2.77], 
[1.97, 2.04], [1.73, 1.79] and [1.72, 1.75] for freshmen, sopho- 
mores, juniors and seniors, respectively. The RO estimated value 
from model was higher than that from epidemiological investiga- 
tion because the close contacts cannot be completely traced. The 
average goodness of fit between the simulation and observation 
data was 76.6% (Appendix figure 1). Also, the estimated Rq value 
for freshmen was higher than those for the other grades. 

Effectiveness of Different Interventions. By using the transmission 
dynamic model, it was estimated that a total of 2,035 cases would 
have resulted by the 20^^^ day if no control measures were imple- 
mented if the outbreak was initiated by a freshman. The total 
number of cases estimated in the same day of the outbreak was 
greater when the outbreak was initiated by a freshman than if the 
outbreak was initiated by a student in any other grade (Figure 4). 

The total number of cases infected by the 20th day of the outbreak 
under different quarantine and isolation situations was estimated 
using the constructed model (Table). Assuming that this outbreak 
was initiated by a freshman, the estimated number of freshmen cases 
under quarantine and isolation was less than the number under the 
quarantine and isolation of equal number of other grades, regardless 
of when quarantine and isolation occurred during the first to fifth 
days of the outbreak (Table; Figure 5). Assuming that this outbreak 
was initiated by a sophomore, the estimated number of cases under 
quarantine and isolation of sophomores was less than the number 
under the quarantine and isolation of equal number of other grades if 
quarantine and isolation occurred on the first, second and third day 
of the outbreak. However, the estimated number of cases under 
quarantine and isolation of freshmen was less than the number under 
quarantine and isolation of equal number of other grades, when the 
quarantine and isolation occurred on the fourth or fifth day of the 
outbreak (Table; Appendix figure 2). Assuming that this outbreak 
was initiated by a junior, the estimated number of cases under quar- 
antine and isolation of juniors was less than the number under quar- 
antine and isolation of equal number of other grades, when the 
quarantine and isolation were performed on the first, second and 
third day of the outbreak. However the estimated number of cases 
under quarantine and isolation of freshmen was less than the number 
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Figure 1 | Number of confirmed HlNl cases in four grades based on the illness onset. 
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Figure 2 | Disease transmission among confirmed cases, (a) Disease transmission between confirmed cases. Red ellipse: freshmen, yellow box: 
sophomores, green diamond: juniors, blue triangle: seniors. The numbers means illness onset, 1 refers to August 27, 2 refers to August 28, and so forth, 
(b) The disease transmission across the grades. I: freshmen, II: sophomores. III: juniors, and IV: seniors. 

under quarantine and isolation of equal number of other grades if was initiated by a senior, the estimated number of cases under quar- 
quarantine and isolation occurred on the fourth or fifth day of the antine and isolation of seniors was less than the number under quar- 
outbreak (Table; Appendix figure 3). Assuming that this outbreak antine and isolation of equal number of other grades, when the 




Grade 



Figure 3 | Comparison of actual Rq values for different grades, t single factor analysis of variance (P < 0.01), the median of actual Rq values for freshmen 
is 2, and is 1 for other grades. 
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Figure 4 | Estimated numbers of cases of A(HlNl)pdm09 when no control measures are used to control the outbreak, varying the source of the index 
case. The model shows predicted cases if the outbreak was initiated by a freshman, sophomore, junior, or senior, respectively. 



quarantine and isolation were performed on the first, second and 
third day of the outbreak. However the estimated number of cases 
under quarantine and isolation of freshmen was less than the number 
under quarantine and isolation of equal number of other grades if 
quarantine and isolation occurred on the fourth or fifth day of the 
outbreak (Table; Appendix figure 4). 

Discussion 

During the characterization of disease transmission among students 
in a college involved in an A(HlNl)pdm09 influenza virus outbreak, 
the data showed that most (93.7%) of the disease transmission chains 
constructed on the basis of laboratory-confirmed cases occurred 
within the same school grade. Cauchemez et al also found similar 
transmission dynamics when investigating an elementary school 
outbreak of A(HlNl)pdm09 influenza among a semirural popu- 
lation in the United States using Bayesian statistical methods^. This 
transmission pattern may be caused by participation in more clus- 
tering activities among students of the same grade. Other school 
outbreak investigations showed that A(HlNl)pdm09 influenza 
spread in the school was likely facilitated by increased student inter- 
actions in the activities, including religious and club activities, close 
contact with friends, and lectures or laboratory work^°'^^'^^. Thus 
interventions to reduce transmission in collective activity may play 
a role in curbing and mitigating future school outbreaks. 

Investigations at the country- or region- level indicated that sea- 
sonal changes in patterns of social contacts have a marked influence 
on the spread of A(HlNl)pdm09 influenza^'^^'^^'^^. In particular, the 
patterns of school terms and holidays affect the incidence^^. Our 
finding of star-shaped rather than branch- and line-shaped trans- 
mission trees with a freshman student as the infectious source 
showed that, even in a micro -environment like this school, patterns 
of social mixing also impacted disease transmission. Freshmen, as 
newcomers to the school and lacking stable social circles, are more 
likely to make new contacts by increasing their peer-to-peer activities 
thus creating their own social networks. The calculated Rq also 
indicated that freshmen facilitated the spread of A(HlNl)pdm09 
influenza between the students. Therefore, infectious disease trans- 
mission modes in the school were established by social contacts of 
freshmen with other students. 

Furthermore, the Rq value for freshmen calculated from the invest- 
igation data or estimated from a model, was greater than those for 
other grades. These findings indicate that freshmen are an important 
determinant of social mixing among students and disease transmis- 
sion was increased by the degree of social mixing between freshmen 
and other students. This was further supported by the model estima- 
tion that many more cases were involved if the outbreak was initiated 
by a freshman than if the outbreak was initiated by a member of 



another grade. Studies using mathematical models also showed that 
patterns of social mixing played a key role in shaping disease 
transmission^'^. 

Studies have shown that many outbreaks of A(HlNl)pdm09 
influenza occurred in schools and school- aged individuals were more 
likely to be infected^°'^^'^^'^«'^«'^^'^^'^^'^^'^^'^^'^^ School-aged individuals 
facilitated the introduction and spread of A(HlNl)pdm09 influenza 
in households''^. It was inferred that the butterfly effect of disease 
transmission in A(HlNl)pdm09 influenza would be observed as the 
disease transmission pattern in school outbreaks; this is influenced 
by social mixing patterns of newcomers in schools. 

During the early stage of A(HlNl)pdm09 influenza pandemic, 
fengxiao (campus quarantine), a reversal of the usual strategy of 
school closures adopted in many other countries, is usually imple- 
mented in mainland China to prevent movement of infected cases 
from the community. This measure was put in place to prevent 
explosive school outbreaks and decelerate disease spread. Studies 
have shown that the early implementation of fengxiao can delay 
the epidemic peak signiflcantly and prevent the disease spread to 
the general population'^'^^. Once a school outbreak has occurred, 
further quarantine and isolation of students is necessary to decrease 
the effect of outbreak without school closure. A key question that 
remains to be answered is how to implement the quarantine and 
isolation measures to effectively mitigate school outbreaks, especially 
in a situation when not all of the contacts can be traced due to 
complex mixing of students in many school activities. 

By using a transmission dynamic model to estimate the case num- 
ber involved in this outbreak under different quarantine and isola- 
tion methods, we found that the quarantine and isolation of 
freshmen more effectively mitigated the school outbreak, if the out- 
break was initiated by a freshman and no matter when quarantine 
and isolation measures were conducted. The same model was used to 
estimate the number of cases if the outbreak was initiated by other 
students in other grades and if quarantine and isolation was con- 
ducted on days two through flve of the outbreak. The model demon- 
strated that quarantine and isolation of freshmen without school 
closure can effectively contain a school outbreak when not all of 
the contacts can be traced. As mentioned above, school- aged indivi- 
duals facilitated the introduction and spread of A(HlNl)pdm09 
influenza in households''^. Therefore, the quarantine and isolation 
of freshmen without school closure is predicted to effectively reduce 
the spread and impact of pandemic influenza in the community. 
From another perspective, this finding suggested that freshmen were 
the key determinant shaping the disease transmission in this school 
outbreak. 

In this study we ignored the heterogeneity of innate immunity to 
A(HlNl)pdm09 influenza among the students, since most people 
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Table 1 | Estimated number of cases involved on the 20th day of this outbreak under the quarantine and isolation of different proportional 
students from different grades on the first day, second day, third day, fourth day and fifth day of the outbreak, assuming the outbreak was 
initiated by different student grades respectively 



First case 



Population intervened 


Proportion 
intervened 




Estimated number of 


cases 






Onset date of 
intervention: 


First day 


Second day 


Third day 


Fourth day 


Fifth day 


Freshmen 


25% 




845 


981 


1025 


1090 


1 136 


Sophomore 


25% 




1 903 


1916 


1931 


1942 


1992 


Juniors 


25% 




1 y4z 


\ y4o 


1959 


1970 


2017 


Seniors 


25% 




1 966 


1 971 


1976 


1982 


2026 


Freshmen 


50% 




609 


~7 A A 
/44 


791 


868 


922 


Sophomore 


50% 




1 842 


1 856 
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1882 


1929 


Juniors 


50% 




1 90o 


1 yUV 
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1981 


Seniors 


50% 




1 929 


1 939 


1949 


1958 


2001 


Freshmen 


75% 




o z z 

o66 


/(OO 
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701 


Sophomore 


75% 




1 776 


1 790 


1804 


1815 


1861 


Juniors 


75% 




1 oo 1 


1 O X. o 
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1894 


1943 


Seniors 


75% 




1 y 1 U 


1 oo 1 
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1932 


1976 


Freshmen 


25% 
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101 1 


1045 


1 109 
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Sophomore 


25% 
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/Ay 


815 
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1 Z 1 1 


1 O 1 o 
1 Z 1 O 
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1280 
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Seniors 
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1 00~7 
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1250 


1282 


1283 


Freshmen 


50% 




880 


899 


927 
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1020 


Sophomore 


50% 




469 


680 


751 


860 


988 


Juniors 


50% 




1 zUU 


1 O 1 "7 
1 Z 1 / 


1221 


1278 
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Seniors 


50% 




1 220 


1 221 


1247 


1280 


1282 


Freshmen 


75% 




/6J 


777 


796 


833 


846 


Sophomore 
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388 


612 


687 


805 


858 


Juniors 
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1191 


1215 


1218 


1276 


1279 


Seniors 
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1214 


1 220 


1244 


1278 


1280 


Freshmen 
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818 


869 


912 


957 


982 


Sophomore 
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101 8 


1 025 


1060 


1 102 


1 107 


Juniors 
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488 


71 8 


786 


856 


1 101 


Seniors 
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1 075 


1 076 


1089 
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1 129 
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1003 


1005 
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1032 


1046 
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Seniors 
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483 


636 


771 


874 


981 


Freshmen 


50% 




715 
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752 


789 
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Sophomore 


50% 




961 


981 


988 


1009 
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Juniors 


50% 
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1029 


1031 


1044 


1100 


Seniors 
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595 


743 
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861 


Freshmen 
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Sophomore 


75% 




937 
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Juniors 


75% 




1015 


1028 
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1040 
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Seniors 


75% 




391 


556 


718 


780 


837 



Freshman 



Sopho 



Junior 



Senior 



were expected to lack innate immunity to the novel strain of influ- 
enza. Additionally, the dynamic model assumed homogenous mix- 
ing of the students from the same grade. Although these limitations 
mean that the mathematical model used in this study represents a 
simplified version of a complex infectious disease system, they may 
provide insight into transmission characteristics of A(HlNl)pdm09 
influenza among students in school outbreaks and serve as a tool for 
decision making when there is uncertainty. We conclude that, dif- 
ferent grades have different transmission modes, and freshmen have 
higher activity and are more likely to cause infection. Freshmen were 
the main facilitators of the spread of influenza during this school 
outbreak, so control measures (e.g. close contact isolation) priority 
used for freshmen would likely have effectively reduced spread of 



influenza in school settings. It will provide the useful suggestions on 
making effective and economic control measures in future school 
outbreak of influenza or other respiratory infectious diseases. 

Methods 

Ethics Statement. The institutional review board of the Academy of Military Medical 
Sciences waived the need for written informed consent from the participants. This 
study is approved by the Academy of Military Medical Sciences Review Board. The 
study is entirely based on epidemiological survey data which is required by Chinese 
Ministry of Health to record during the epidemic control. There was no request for an 
ethical permission according to Chinese law. Even so, we took the following steps to 
protect our survey respondents and anyone mentioned in the survey. Firstly, we kept 
all completed questionnaires under strict anonymity and confidentiality, left no 
identifiable personal information (such as names of respondents or anyone listed in 
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Figure 5 | Estimated numbers of cases under quarantine and isolation of equal freshmen, sophomores, juniors, or seniors over the first 20 days of the 
outbreak, if the outbreak was initiated by freshman, (a), quarantine and isolation of 25% freshmen, sophomores, juniors, or seniors on the first day of 
outbreak; (b), quarantine of 50% freshmen, sophomores, juniors, or seniors on the first day of outbreak; (c), quarantine of 75% freshmen, sophomores, 
juniors, or seniors on the first day of outbreak; (d), quarantine of 25% freshmen, sophomores, juniors, or seniors on the second day of outbreak; 
(e), quarantine of 50% freshmen, sophomores, juniors, or seniors on the second day of outbreak; (f), quarantine of 75% freshmen, sophomores, juniors, 
or seniors on the second day of outbreak; (g), quarantine of 25% freshmen, sophomores, juniors, or seniors on the third day of outbreak; (h), quarantine 
of 50% freshmen, sophomores, juniors, or seniors on the third day of outbreak; (i), quarantine of 75% freshmen, sophomores, juniors, or seniors on the 
third day of outbreak; (j), quarantine of 25% freshmen, sophomores, juniors, or seniors on the fourth day of outbreak; (k), quarantine of 50% freshmen, 
sophomores, juniors, or seniors on the fourth day of outbreak; (1), quarantine of 75% freshmen, sophomores, juniors, or seniors on the fourth day of 
outbreak; (m), quarantine of 25% freshmen, sophomores, juniors, or seniors on the fifth day of outbreak; (n), quarantine of 50% freshmen, sophomores, 
juniors, or seniors on the fifth day of outbreak; (o), quarantine of 75% freshmen, sophomores, juniors, or seniors on the fifth day of outbreak. 



the contact diaries) in our electronic data set, and destroyed the completed 
questionnaires after data cleaning. Secondly, we used the personal information from 
our sample list only to locate and identif)^ our targeted respondents and then 
destroyed the name list after the interviews. 

Outbreak Survey. The outbreak began on a college grounds near Beijing in northern 
China on August 27, 2009 and an epidemiological investigation was conducted by the 
team from the Institution of Disease Control and Prevention. The college awards a 
four-year degree; all students reside at the college with a total enrollment of 6721 
students (1742 freshman, 1805 sophomores, 1808 juniors, 1366 seniors) ranging in 



age from 19 to 24. Influenza"like illness (ILI) was defined as the sudden onset of a 
fever of >38°C with cough or sore throat in the absence of other diagnoses; a 
confirmed case of influenza was defined as a case of ILI with laboratory" confirmed 
A(HlNl)pdm09 influenza virus infection using the polymerase chain reaction (PCR) 
according to WHO guidelines'^. We detected the antibody levels against seasonal 
influenza A (HlNl), H3N2 and influenza B, and found different grades had no 
difference. Each confirmed case was interviewed face"tO"face to obtain demographic 
information, clinical symptoms, onset date of influenza"like symptoms and contacts. 
A close contact was defined as a person having cared for, lived with, or who had direct 
contact with respiratory secretions or body fluids of a probable or confirmed case, or 
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conversation with probable cases more than 20 minutes continuously in one 
encounter within 3 days after symptom onset^l Public health control measures 
implemented by the response team, led by the Institute of Disease Control and 
Prevention of the Academy of Military Medical Sciences, followed the guidelines from 
the Ministry of Health of China. From September 1, interventions, including 
individual isolation of confirmed cases and quarantining of contacts (home stay or 
dormitory) were taken. All students and staff were not permitted to leave the school 
grounds. Interventions were sustained during the whole period of the epidemic. 
Before the intervention measures were taken, all four grades had confirmed cases. 

Transmission Chains. We constructed disease transmission pathways among the 
confirmed A(HlNl)pdm09 influenza cases. If student A and student B are close 
contacts and symptom onset appeared in student A prior to student B, we assume that 
student A infected student B. We used the expression "A B" to show the 
transmission relationship between A and B. If there were several possible routes of 
transmission, we assumed the "earlier infector" is the "real infector". Disease 
transmission pathways were constructed using the program Pajek^^. Pajek is a 
Windows program for analysis and visualization of large networks having some 
thousands or millions of vertices. 

Transmission Dynamic Model. The susceptible-exposed-infected-removed (SEIR) 
model was used to simulate transmission dynamic processes. The SEIR model has 
been widely used in epidemic simulations The SEIR model was used to simulate 
infectious diseases transmission dynamics, estimate epidemiological parameters and 
effectively evaluate interventions^^"^*^. The SEIR model separates the population of a 
region into four compartments, susceptible, exposed, infected, and removed, and 
simulates the time evolution of each of these subpopulations. We assumed the school 
to be a closed setting and the students from the same grade (e.g. freshmen, 
sophomores, juniors, seniors) to be homogenously mixed in school activities due to 
students from same grade with almost same age, and also they have same schedule, 
and ignored the heterogeneity of innate immunity to A(HlNl)pdm09 influenza 
between the students. Students from each grade were divided into four groups 
(Susceptible, Exposed, Infected and Removed (or Recovered)). Under these 
circumstances, a transmission dynamic model was constructed (Figure 6). 

There were M individuals in grade i, i = 1, 2, 3, 4. Therefore the number of 
susceptible students from grade / in the t day of outbreak; is denoted as Si(t); Ei{t), is 
the number of exposed students from grade / in the day t of outbreak; Ii{t), is the 
number of infected students from grade i in the t day of outbreak; and Ri{t) is the 
number of recovered students from grade i in the t day of outbreak. We also denoted 

(i = 1, 2, 3, 4), as the average number of secondary cases per primary case in a 
susceptible population from the same grade in the infectious period of primary cases 
per unit time; ocg, as the average number of secondary cases per primary case in a 
susceptible population from the different grade in the infectious period of primary 
cases per unit time. A parameter (C) was set to control the change of (i = 1,2,3,4,5) 
under the control measures taken. The parameter C is equal to 1 when no inter- 
ventions are taken, it is smaller than 1 when interventions are taken, and the more 
stringent the control measures, the smaller the C value. From the literature- derived 
parameter values of the incubation and infectious periods\ the incubation period was 
assumed to be 1/p =1.1 days and the infectious period was l/r| = 2.5 days. The ai/r| 
is the RO value of grade i (i = 1,2,3,4). The newly infected cases per day were used to fit 
our model. The ANNEAL algorithm was used to optimize parameters of function, 
which is a stochastic method currently in wide use for difficult optimization pro- 
blems^'*'^^. The model was constructed by using ordinary differential equations^^: 
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Using the constructed model, the initial dynamics of disease transmission among 
the different grades were reproduced, and used to estimate Rq values for freshman, 
sophomores, juniors or seniors We simulated 100 times using the ANNEAL algo- 
rithm that minimizes a function using the method of simulated annealing^^. 
Simulated annealing (SA) is a generic probabilistic metaheuristic for the global 
optimization problem of locating a good approximation to the global optimum of a 
given function in a large search space. We simulate effectiveness of intervention 
through control numbers of four compartments (Susceptible, Exposed, Infected, and 
Removed). 

Estimation of Intervention Effectiveness. By using the constructed model, estimates 
were derived for the total numbers of cases involved on the 20th day of the outbreak. 
Assuming that the first case of this outbreak was a freshman, estimates were made of 
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Infected 



Removed 
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Removed 



Figure 6 | The transmission dynamics model of A(HlNl)pdm09 
influenza among students attending a college in northern China, t i = 

1,2,3)4, j = 1,2,3,4, i j. {Each grade was divided into four groups 
(Susceptible, Exposed, Infected and Removed). 

the number of infected students involved in the outbreak on the 20th day, when 25%, 
50%, or 75% of freshmen, sophomores, juniors, or seniors were quarantined and 
isolated on the first, second, third, fourth and fifth day of the outbreak. Similar 
estimates were made assuming the first case was a sophomore, junior, or senior, 
respectively. For these estimates, we applied the same stringent intervention measures 
to contain the epidemic of A(HlNl)pdm09 influenza, including isolation of 
confirmed cases and quarantining of contacts, which were recommended by the 
Chinese Ministry of Health on May 29, 2009^. The susceptible students involved in 
this outbreak were assumed to be close contacts. Interventions were sustained 
through the whole epidemic in order to avoid creating a population of individuals 
who were protected from infection but would become infected after the intervention 
was stopped. The effectiveness of different quarantine methods under the different 
conditions were compared according to the estimated number of cases involved. 
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